Mon. Not. R. Astron. Soc. 000, [THl (2009) Printed 29 August 2009 (MN lAT^ style file v2.2) 



Infrared identification of 4U1323-619 revisited 



Ivan Yu. Zolotukhin^*, Mikhail G. Revnivtsev^'^, and Nikolai I. Shakura 

^ Sternberg Astronomical Institute, Moscow State University, Universitetskij pr., 13, 119992, Moscow, Russia 
'^Excellence Cluster Universe, Technische Universitdt Miinchen, Boltzmannstr.2, 85748 Garching, Germany 
Space Research Institute, Russian Academy of Sciences, Profsoyuznaya 84/32, 117997 Moscow, Russia 



Accepted 2009 August 29. Received 2009 August 29; in original form 2009 August 11 



ABSTRACT 

We re-examine the infrared counterpart of the dipping low-mass x-ray binary 4U1323- 
619. New X-ray data available from the XMM and Chandra observatories combined 
with archival IR observations from the ESO 3.6m New Technology Telescope allow us 
to define a new possible counterpart. We present here its photometric properties and 
compare them with a simple analytical model of an accretion disc illuminated by the 
hot central corona known to be present in the binary system. 
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1 INTRODUCTION 

■ 4U1323-619 is a dippi ng low-mass x-ray binary detected 
' by U huru and Ariel V (jPorman et al.l[l978l : IWarwick et al.l 
I Il98ll ) . It exhibits X-ray bursts and irregular intensity dips 

■ repeating every 2.94 h. Dips discovered in EXOSAT data 

■ (jParmar et al.l ll989l ) lasting typically 30 per cent of the 
, orbital cycle, and the lack of X-ray eclipses are indicative 

■ that the binary pla ne is viewed at inclination i = 60 — 80° 
l|Frank et al.l Il987h . Such dipping behaviour is explained 

. by periodic obscurations of the central X-ray emitting re- 
' gion by a structure locat ed in the outer regions of the 
; disc (| White fc Swankll982h . High interstellar matter column 

■ density (Nh ^ 4 x 10^'^ cm ~'^) derived from X-ray spectral 
' models (jParmar et al.l[T989h suggests that this source pos- 
sesses significant extinction of the order of Ay ~ 1 8 mag 
and, t hus, cannot be observed at optical wavelengths. ISmald 
l|l995l ) attempted to find its IR counterpart subject to much 
lower extinction and reported a candidate which demon- 
strated some variability of the IR flux (though at a low 
significance level). 

Similarly to other sources, the dipping activity 
in 4U1323-619 attracted attention of researchers as 

a laborator y for study ing the a ccretion disc corona 

(iBalu cihska- Church et alJ Il999l : iBoirin et all l2005l : 
IChurch et all |2005| ). These studies and other recent 
observational efforts produced new data available for this 
source in archives of X-ray missions and the European 
Southern Observatory (ESO). We have noticed a significant 
positional discrepancy between the Einstem coordinates of 
4U1323-619 and its more recent XMM and Chandra obser- 
vations. Having found the IR archival data for 4U1323-619 
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by the means of the Virtual Observator y, we decided to 
re-examine the IR candidate discovered bv lSmald (|l995l ). 



2 OBSERVATIONS AND RESULTS 
2.1 X-ray observations 

We have analysed Chandra observations performed on Sep 
25, 2003 with a total exposure time about 40 ksec. The ob- 
servations were aimed at high energy resolution spectroscopy 
and the telescope was equipped with gratings. Detectors 
were operating in the so-called Continuous-Clocking (CC) 
mode, when the information in one spatial direction is lost. 

Therefore we have determined the source coordinates 
mainly in one direction (with the nominal Chandra lo- 
calisation accuracy ~ 0.6 arcsec at 90 per cent confi- 
dence). Orientation of the effective position "strip" de- 
pends on the orientation of the telescope during the ob- 
servations. We represent the Chandra localisation by a 
wide ellipse in Fig. [T] The centre of the ellipse is at 
RA(J2000)=201.65397 deg, Dec(J2000)=-62.135403 deg, its 
positional angle is 156.6 deg, the ellipse is ~ 0.6 arcsec wide 
in one direction and effectively infinitely long in another. In 
Fig. [T] we have limited the size by arbitrarily adopting a 
length of 5 arcsec. 

XMM-Newton observed the source position two times, 
on Aug 17, 2001 and Jan 29, 2003, the observations span- 
ning 26 ksec and 51 ksec, respectively. Data of the XMM- 
Newton EPIC-MOS cameras were analysed using the stan- 
dard tasks of the Science Analysis Software (SAS) vQ.O.O. 
Unfortunately, the source is sufficiently bright to have a non- 
negligible effect on the XMM/EPIC-MOS imaging capabil- 
ity because of pile-up. This depresses the brightness of cen- 
tral pixels, degrading the localisation accuracy, in spite of 
the possibility of cross-calibrating the absolute astrometry 
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Table 1. Positions and Ks magnitudes of the sources from Fig.[T] 
The astrometric uncertainty is 0.2 arcsec for all objects. Source 
B is the proposed new IR, counterpart for 4U1323-619, marked 
boldface. 



Source 


RA (J2000) 


Dec (J2000) 


Ks 


A 


13:26:37.22 


-62 08 07.7 


18.56 ±0.23 


B 


13:26:36.98 


-62:08:09.0 


18.12 ±0.20 


C 


13:26:36.25 


-62:08:10.1 


17.87 ±0.18 



magnitudes and the astrometric positions which are given 
in Tab. [1] The single object within the Ei nstein error cir- 
cle (source C in this work) was proposed bv lSmal3 l| 19951 ) as 
the IR counterpart of 4U1323-619 and can now be ruled out. 
This mis-identification was probably caused by the under- 
estimation of the position uncertainty radius for the source, 
which is quite faint for Einstein /¥LKl. We propose the only 
source visible inside Chandra and XMM positional errors in- 
tersection down to a limiting magnitude Ks ~ 19.4 mag (So- 
upper limit), source B, as a new counterpart for 4U1323-619. 



of XMM using optical sources within its field of view. There- 
fore, we would like to adopt some conservative uncertainty 
radius of the source position ~ 3 arcsec approximately cor- 
responding to 90 per cent confidence, which without pile- 
up might have been significantly improved. The source po- 
sition was determined to be RA(J2000)=201. 654870 deg, 
Dec(J2000)=-62.135985 deg. 



2.2 IR observations 

The infrared data analysed in this work were obtained on 
May 18, 2005 between 03:40 and 04:01 U T using the SOFI 
infra red spectrometer and imaging camera (jMoorwood et all 
1 19981 ') at the ESO 3.6m NTT telescope under programme 
ID 075.D-0529(A). All scientific and calibration data were 
retrieved through the publicly available ESO observational 
archivcQ. The scientific data contain 21 three second expo- 
sures, each taken in Small Field imaging mode with pixel 
scale 0.144 arcsec pix~^ and seeing ~1.1 arcsec. Due to 
the specific nature of IR observations, i.e. a rapidly chang- 
ing background and the techniques used to determine this, 
the reduction procedure is quite different from optical imag- 
ing observations. To get the data ready for scientific analy- 
sis we, therefore, used tools and recipes provided by ESO: 
GASGANO v2.3.0ior data organisation tasks and the SOFI 
data reduction recipes from the Common Pipeline Library 
v5. 0.0 to correct for bias, flat fleld and frame jittering. After 
co-adding 21 reduced co-aligned frames (see Fig.[l}, we per- 
formed aperture photometry measureme nts with the SEx- 
tractor software (|Bertin fc Arnoutil 19961 ) and calibrated in- 
strumental magnitudes using 2MASS Ks photometry of field 
stars. Our photometric errors therefore include calibration 
dependency uncertainties of 0.10 mag. The astrometric so- 
lution was obtaine d in the 2MA SS reference frame with the 
SCAMP software (iBertinlbOO^ ) and had 0.2 arcsec calibra- 
tion uncertainty. Infrared measurements of all sources of in- 
terest in the 4U1323-619 field are listed in Tab.[Tl 



2.3 Results 

We overplotted positional uncertainties obtained from the 
XMM and Chandra data on the combined IR image (see 
Fig -Ul, to g ether with the Einstein error circle which, follow- 
ing 'Smalc (1995), is centred at the X-ray position of Source 
D from Parmar ct al. (1989), as measured by Emstetn HRI: 
aj2ooo = 13'^26'"36.08% Jjaooo = -62°08'10.2", R = 
2.5 arcsec (90 per cent confidence). In Fig. [1] we designated 
sources of the interest by a capital letters and measured their 
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3 DISCUSSION 

We obtained simple estimations of optical luminosity of the 
accretion disc illuminated by the central isotropic spheri- 
cal X-ray sou rce with a 0.1-10 k eV luminosity Lx ~ 5.2 x 
10 erg s' iBoirin et al.l 120051 ) assuming for convenience 
that its radius is ~ few x lORg which means that it is rel- 
atively small comparing to the accretion disc size. While it 
is known th at the source has varied in luminosity over a 20 
year period (|Baluciiiska-Church et all[200& ) and an interpo- 
lation to the date of the IR observations gives ~ 2.5 times 
smaller luminosity, we leave mentioned Lx value because 
the bolometric luminosity is expected to be the same factor 
higher. Given a system period P = 2.94 h and assuming a 
compact object to be a neutron star with Mx = IAMq, one 
can estimate the mass of the secondary component from the 
mass-radius relation for main sequence stars and the con- 
strain t that the star fills its Roche lobe of size (jEggletonI 
Il983l) : 



r/a = 0.49 



^2/3 



0.6<j2/3 +ln(l + gl/3) 



(1) 



where q — Mopt/Mx and a is the semi-major axis of the 
binary. This simple estimate gives Mopt — 0.25 M(7) and en- 
ables us to determine the accretion disc size using IPaczvnskil 
(1977): rout ^ 0.51a = 0.62i?Q for fj, = Mopt/{Mopt + Mx) = 
0.15. 

The effective temperature of an accretion disc with a 
height zo cx r" ill uminated by a point source in its centre is 
approximated by (|Shakura fc Sunvaevll 197*3 ): 



Lx „ Lx I T ^ 20 



^, (2) 
47rr^ VKr^ r 

where r] is the absorbed fraction of the radiation impinging 
on the disc surface assumed to be ~ 0.5. Let us consider the 
situation when the illumination does not change much the 
structure of a standard a— disc. We therefore adopt stan- 
dard Q-disc parameters, n = 9/S:,ZQ/rout — 0.05, giving a 
temperature Te//, out = 10500 K at the outer disc radius 
r'out ■ 

The observed flux (from a single side of the accretion 
disc as visible from Earth) at frequency 1/ = 1.4 X 10" Hz 
for the Ks fllter is close to the Ray leigh- Jeans range and can 
be estimated as follows: 



^ 64 . 2v'kT,ff, outrLt - 
if = — TTCOsi „ ,„ , e 



49 



(3) 



-15/32 



Here the factor 64/49 accounts for the Teff{r) cx r 
dependency in the standard accretion disc, and / ~ 1.3 
compensates for a non Ray leigh- Jeans law for a given 
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Figure 1. The co-added Ka image of ttie 4U1323-619 field constructed from 21 individual three second exposures after reduction and 
correction for frame jittering. Einstein, Chandra, and XMM positional uncertainties are overplotted. North is up, East is left. Measured 
sources are designated by capital letters; see the text for details. 



hu/kTsff, out ~ 0.58 ratio. Moreover, the observed flux is 
obviously decreased by the interstellar extinction with op- 
tical depth T^. Taking a moderate orbital plane inclination 
angle i — 70°, a source distance d = 10 kpc, and a typi- 
cal extinction Ak,, — 2.0 mag (r^ ~ 1.84) derived from the 
3D g alactic extinction map in this direction (,Marshall et al.l 
l200d) . gives = 2.1 X 10"^^ erg cm"^ s"^ Hz"^ which is an 
order of magnitude less than the observed flux (3.1 .. . 4.5) x 
10~^* erg cm~^ s~^ Hz~^ translated from the measured 
magnitude Ks = 18.12 ± 0.20 mag a nd its uncertainty using 
zeropoints from I Cohen et al] l|2003l '). 

We also estimated the observed flux in another limiting 
case for an illuminated disc with constant temperature along 
the z coordinate (i s other mal disc). Following equation (A5) 
from IVrtflek et al] l|l990l ) with the same set of the system 
parameters used above, one gets Te//, out = 10000 K. Use of 
Teff(r) oc r~^'^^ in integrating across the disc for isothermal 
model gives a factor of 14/11 instead of 64/49 in Eq. |3]and 
the observed flux becomes a factor of 0.9 smaller than 
in the flrst case, considering also slightly different Te//, out- 

There are several possible explanations for this discrep- 
ancy between the observed and predicted fluxes. The ex- 
istence of a hot (2 X 10^ K) atmosphere ab ove the outer 
regions of a disc (jjimenez- Garate et al]|2002f ) can mcrease 
the observed flux at a given frequency u = lA x 10^'* Hz 
by a factor of 1.3 . . . 1.7 if we naively consider the effect of 



increased zo/rout — 0.2 ratio on Eq. [2H31 In such a hot at- 
mosphere, the scattering of X-ray photons on free electrons 
with consequent penetration to sub-photosphere layers and 
thermalization takes place. This increases disc illumination 
and hence brightens the optical emission. 

Also we assumed in the beginning the compact na- 
ture of the central X-ray source while it is now known that 
there exists extended accretio n disc corona with the radius 
300 00 km for adopted L x (IChurch fc Balucihska-ChurchI 



'2004'; 'Bal uciriska-Church et al.ll2009l ). This changes the il- 
lumination geometry of the outer regions of a disc, namely 
cos6' ~ zo/r in Eq. [2j and thus increases Te//, out and Qi, 
by a factor of (n — ~ 1.7, correspondingly. 

Placing the obj ect to the distance o f 4-5 kpc instead of 
10 kpc suggested bv lParmar et al.l l| 19891 ) on the basis of ob- 
served bursts being sub-Eddington would account for all in- 
consistency of predicted vs. observed fluxes because of signif- 
icantly reduced extinction in this direction and the fact that 
Qi, scales with distance as dT^^'^ for a fixed X-ray flux (which 
follows immediately from Eq.[2H3]since oc Te//, out/<f oc 
iJ^Jr^ jd} while itself Lx oc for the mentioned condition). 



2.0x10 



'^^ erg cm s 



^Hz- 



We therefore might get 
in this case. 

There is also a possibility of existence of a ~ 2 — 3 mag 
fainter object inside the combined XMM and Chandra error 
box which could be the actual counterpart of 4U1323-619, 
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but this cannot be ruled out on the basis of existing archival 
data and requires dedicated high spatial resolution X-ray 
observations and deeper phase-resolved follow-up IR obser- 
vations, probably with adaptive optics in order to reduce 
field contamination by the source B. Discovery of the 2.94 h 
variability period in the IR source would allow to finally 
identify the 4U1323-619 counterpart. 



4 CONCLUSIONS 

On the basis of examination of archival ESO NTT data 
within the area of overlap of XMM and Chandra error re- 
gions, we have identified a probable candidate for the IR 
counterpart of 4U1323-619. Its observed Ks magnitude sig- 
nificantly differs from that predicted by a simple analyti- 
cal model of an accretion disc illuminated by a hot central 
spherical corona with a parameter set available in the lit- 
erature. While putting the object to 4-5 kpc instead of as- 
sumed 10 kpc would explain the discrepancy, wc encourage 
high spatial accuracy X-ray observations (e.g. using Chan- 
dra/HRC) and deeper phase-resolved follow-up observations 
inside the resulting error box listed here. 
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